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1. Abstract 


A review of the important constraints on gravity wave induced diffusion 
of chemical tracers, heat, and momentum is given. Ground-based microwave 
spectroscopy measurements of H 2 0 and CO and rocket-based mass spectrometer 
measurements of Ar constrain the eddy diffusion coefficient for constituent 
transport (K zz ) to be (l-3)x 10 5 cm 2 s* 1 in the upper mesosphere. Atomic 
oxygen data also limits !(„ to a comparable value at the mesopause. From the 
energy balance of the upper mesosphere the eddy diffusion coefficient for heat 
transport (D H ) is, at most, 6 x 10 5 cm 2 s' 1 at the mesopause and decreasing 
substantially with decreasing altitude. The available evidence for mean wind 
deceleration and the corresponding eddy diffusion coefficient for momentum 
stresses (D„) suggests that it is at least 1 x 10 6 cm 2 s _1 - in the upper 

mesosphere. Consequently the eddy Prandtl number for macroscopic scale 
lengths is > 3 . 
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2. Introduction 


The importance of gravity waves in the mesosphere and lower thermosphere 
has been known for more than twenty five years. The seminal paper by Hines 
(1960) opened up an area of research which unfortunately took over twenty 
years to realize the full significance of gravity waves in this region of the 
atmosphere. Research efforts by Pitteway and Hines (1963), Lindzen (1967, 

1968), and Hodges (1969) provided the needed foundation for present day 
research on gravity wave induced diffusion, but did not then yield a 
parameterization of the gravity wave stresses on the mean circulation and 
structure of the middle atmosphere. It was not until 1981 when Lindzen (1981) 
published his perceptive recipe for parameterizing the gravity wave momen 
stresses on the mean zonal flow that the full impact of gravity waves on the 
middle atmosphere was comprehended. An unfortunate result of this important 
paper was the almost sole emphasis on the momentum balance of the middle 

atmosphere . 

Dynamical models were constructed which examined in great detail the 
deceleration effects of gravity waves on the zonal flow but omitted any 
thermodynamic effects on the global mean temperature structure (e.g. Holton, 
1983; Garcia and Solomon, 1985). However in the Garcia and Solomon (1985) 
study diffusive effects on temperature departures from the global mean 
temperature were included. The consensus of these studies was that the eddy 
diffusion coefficient for momentum stresses (Dm) is - 10 cm s (Lindzen, 

1981; Holton, 1982, 1983; Garcia and Solomon, 1985). It should be noted that 
this quantity is proportional to , where a is the mean zonal wind and c 

is the phase speed of the gravity wave. This extreme sensitivity of D M to (u-c) 
makes an accurate calculation of its value very difficult from observed values 
of S and c. Because the mean circulation and associated transport of 
constituents and heat are also directly dependent on D u /(U-c) , theoretical 
prediction of their strengths is also subject to the uncertain magnitude of 
(5-c) 1 through the uncertainities associated with the computed radiative drive 
and the input gravity wave spectrum of phase speeds. 
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Johnson and Wilkins (1965) noted that the observed lower thermospheric 
temperature gradient was inconsistent with molecular conduction of heat alone 
and concluded that eddy transport of potential temperature was required. In a 
more quantitative study Johnson and Gottlieb (1970) inferred from the globally 
averaged heat balance that the eddy diffusion coefficient for heat transport 
(D h ) was 1 x 10 5 cm 2 s' 1 at 60 km and increased to 1 x 10 7 cm 2 s' 1 at 120 km. 

Of more importance to us are the values in the 75-85 km region which they 
deduced to be (4-7) x 10 5 cm 2 s' 1 . Hunten (1974) argued that these analyses 
ignored an essential aspect of the physics, namely that dissipation of 
mechanical energy to generate turbulence might yield net heating rather than 
net cooling of the mesopause region. His inference of net heating depended 
crucially on the assumption that turbulence persists only if the Richardson 
number of the flow (Ri) remains at its critical, onset value, generally s 0.2. 
Johnson (1975) noted that turbulence can persist in a flow with Ri as large as 
1 and thus argued that steady- state turbulence would maintain the very stable 
lower thermosphere near Ri — 1, which implies net cooling by the action of 
turbulence on the mean circulation. 

In a similar vein, chemical tracers bear signatures of transport effects 
and Colegrove et al. (1965, 1966) used the 0/0 2 density ratio to infer the 
vertical eddy diffusion coefficient for tracer transport (K rx ) . Their average 
value in the 80-120 km region was in the range of (0.8-8) x 10 6 cm 2 s' 1 , with a 
preferred value of 4 x 10 6 cm 2 s' 1 . It is standard practice in ID 
photochemical models with eddy and molecular diffusive transport to 
empirically deduce the eddy diffusion coefficient profile with a minor 
constituent whose density profile is extremely sensitive to the adopted 
values. In the one extreme, Hunten and Strobel (1974) used argon measurements 
by von Zahn (1970) to deduce a homopause value of ^ - 3 x 10 5 cm 2 s" 1 . In 
the other extreme Keneshea and Zimmerman (1970) argued for highly structured 
profiles of K zz with peak values in excess of 1 x 10 7 cm 2 s' 1 . In the lower 
thermosphere the uncertainty in the heating and cooling rates at that time 
could not exclude such large values of the eddy diffusion. In a turbulent 
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atmosphere, tracer and potential temperature transport by eddies should be 
described by the same coefficient, i.e., D H - K„ , in the limit of a 
chemically inert species and adiabatic motion. For chemically active species 
one must include chemical acceleration of vertical transport which arises 
physically because the mixing ratio perturbation for chemical species with 


finite chemical loss undergoes a phase shift relative to the wave velocity 
fields such that the eddy fluxes are nonzero rather than in quadrature, as in 
the limit of no chemical loss (Strobel, 1981). Strobel et al. (1987) in a 
study of mesospheric chemistry found this term to be potentially important for 

only odd oxygen and then in a subtle way. 

The focus of this paper is on gravity wave induced diffusion over 
macroscopic vertical length scales on the order of the atmospheric scale 
height, H. There is a significant body of literature on measurements and 
interpretation of smaller scale diffusion observed by rocket experiments and 
radars and reviewed in depth by Hocking (1985) . Turbulent layers observed by 
these techniques tend to be "sub -macroscopic" . There is the difficult 
question on how to relate diffusion on these smaller scales to the larger 
scales of interest here. The reader should keep in mind that additional 
constraints may be supplied by this broad data base, although the 
extrapolation to larger scales is not obvious and straightforward. 


3. Constraints on K„ 

The momentum stresses of saturated gravity waves are largest in the upper 
mesosphere and it is there that independent constraints on gravity wave 
diffusion are most important. Investigations by Allen et al. (1981) and 
Strobel et al. (1987) identified H 2 0, CO, 0 3 , and 0 as the chemical species 
most diagnostic of tracer transport in the upper mesosphere. Of these species 
Strobel et al. (1987) argued that the vertical H 2 0 profile was best suited to 
infer the correct magnitude of vertical mixing in the mesosphere. 



In their study they derived an approximate expression for the vertical 
eddy diffusion coefficient in terms of the eddy diffusion coefficient for heat 
transport and chemical acceleration of vertical transport 

^«(0 = D h + 2x 1 0 1 0 Z. , yl (1) 

with the numerical value in cgs units, A is a measure of gravity wave 
amplitude and equal to 1 at the saturation limit, and L x is the chemical loss 
rate of specie i. According to recent measurements by Vincent and Fritts 
(1987) , gravity wave amplitudes vary throughout the year in the range of A - 
0.6-2. With the exception of odd oxygen, all other species which undergo 
significant transport in the mesosphere have L £ A « 5 x 10' 9 D H and thus K zz (i) 

“ d h . 

In Fig. 1 theoretical H z 0 mixing ratio profiles from an ID eddy 
diffusion, photochemical model (Allen et al., 1981 with updates discussed in 
Strobel et al., 1987) are compared with observed values. The corresponding 
vertical eddy diffusion coefficient (D H ) profiles are illustrated in Fig. 2. 

The mesospheric water vapor measurements were obtained by ground-based 
microwave spectroscopy from the Jet Propulsion Laboratory in Pasadena, 
California during the time period April to June 1984 (Bevilacqua et al . , 

1983). At this time vertical motions associated with the mean meridional 
circulation would be expected to be upward (Garcia and Solomon, 1985) and the 
inferred values of D H (-K zz ) from the observed H 2 0 mixing ratio profiles should 
be upper limits. The water vapor is most sensitive to transport in the 70-80 
km region and there the observations clearly constrain D H to a range of (1-2) 
x 10 5 cm 2 s 1 as Models B, C, and D show in comparison to Model A. Model A was 
preferred by Allen et al. (1981, 1984) in order to satisfy the water vapor 
measurement by H. Trinks (private communication, 1979) at 90 km. Ground-based 
microwave spectroscopy measurements imply that the H 2 0 mixing ratio at 80 km 
is comparable to the Trinks' value at 90 km. The rapid dissociative loss of 
H 2 0 in the 80-90 km region requires a basic incompatibility between these 
measurements. If one examines the microwave signature of H 2 0, the signal to 
noise increases with increasing altitude with the highest quality data in the 
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70-80 tan region (Bevilacqua et al. , 1987a). The associated error bars clearly 
exclude the observation of a water vapor profile even remotely close in mixing 

ratio to Model A in Fig. 1. 

Ground-based microwave spectroscopy measurements of CO (Bevilacqua et 
al., 1985) confirm the deduction of slow vertical mixing in the mesosphere 
(Strobel et al. , 1987). But it must be kept in mind that the large vertical 
scale height of the CO density profile makes it somewhat insensitive to 
diffusive transport and hence more susceptible to advective transport. This 
follows from a comparison of the respective time constants (H 2 /K„) and (H/ w) , 
where 5 is the zonally averaged vertical wind. In addition the Solar 
Mesosphere Explorer (SME) near- infrared spectrometer measurements of absolute 
0 3 concentrations are best understood with slow vertical mixing, which yields 
low H 2 0 mixing ratios and odd hydrogen densities. This produces red 
catalytic destruction of 0 3 by odd hydrogen and hence high ozone mixing ratios 
in better agreement with SME measurements. Brasseur and Offermann (1986) 
analyzed 0 concentration measurements and concluded that the vertical eddy 
diffusion coefficient is about 10’ cm 2 s' 2 at the mesopause, also consistent 
with the above results. Note that 0 unlike CO has a very small scale height 
at and below the mesopause which renders it extremely sensitive to diffusive 
transport. It is interesting to note that the argon measurements that guided 
Hunten and Strobel (1974) to adopt low vertical mixing in the mesopause region 
are in excellent agreement with inferences from other species. Evidence from 
absolute concentrations of chemical tracers of mesospheric transport thus 
suggest that (-D H , for all tracers but 0 and 0 3 ) is - (l-3)xl0 5 cm s in 

the upper mesosphere. 

4. Constraints on D H 

As alluded to in the Introduction it is more difficult to obtain powerful 
constraints on the eddy diffusion coefficient for heat transport because it 
involves extracting a small residual from large terms in approximate balance 
in the thermodynamic heat equation. The importance of diffusive transport of 
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potential temperature by breaking gravity waves was demonstrated by Schoeberl 
et al. (1983) in a numerical study of gravity wave breaking and stress in the 
mesosphere. Chao and Schoeberl (1984) emphasized that it is the turbulence 
created by the breaking gravity wave that transports potential temperature 
rather than coherent, HTT heat transport by the gravity wave. Apruzese et 
al. (1984) in a study of the globally averaged temperature of the mesosphere 
and lower thermosphere concluded that D H must be less than 10 6 cm 2 s' 1 . 

Strobel et al. (1985) reiterated the arguments of Apruzese et al. (1984) and 
argued that the eddy Prandtl number (P r - D„/D H ) over macroscopic scale 
lengths must be large if D„ exceeds 10 6 cm 2 s' 1 , because D,S6xio'cm 1 « ‘. 

To quantify this discussion, let us examine the thermodynamic heat 
equation for the globally averaged temperature field, <T> , with < > denoting 
global average and overbar zonal average 


a<7> 

Jt 


pc 


1 1 

' Q l iv ~ Q ir > + 


P C p d Z 


dT i 

> — \<pwT> + <H -C, 

pc p \dz) ^ 0 « 


( 2 ) 


where Q„y is the solar heating rate due to 0 2 and 0 3 absorption, is 
infrared cooling rate, p is the mass density of the atmosphere, i is the 
molecular heat conductivity, c p is the specific heat at constant pressure, 5 
is the zonal ly averaged vertical velocity, 


N 2 D u 


* 2c 


(H-f; 1 ) 


(3) 


is the conversion rate of wave energy to heat with efficiency c, and 


C ’ = ^M PC ' D ‘{£ + c ")) 

is the cooling rate due to the divergence of the downward turbulent or eddy 
heat flux. Here 


n 2 


T\dz c p ) 


(S) 
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is the buoyancy frequency. Note that when D-.'tD.m is less than H the heat 
flux is divergent and there is cooling, whereas if it is greater than H' the 
heat flux is convergent and heating occurs. The expressions for H g and C g 
were derived by Schoeberl et al . (1983) under the assumption of Lindzen (1981) 

that when gravity waves break or saturate their amplitudes remain constant 

with altitude. 

Equation (2) has been solved for the globally averaged, steady-state 
temperature without the term <p*T> which represents the convergence of the 
downward heat flux associated with the anticorrelation of the zonal mean 
vertical velocity (a result of gravity wave breaking) and the zonally averaged 
temperature. For a Prandtl number of 1. Fig. 3 shows illustrative globally 
averaged temperature profiles with the corresponding vertical eddy diffusion 
coefficient profiles in Fig. 4. These results led Apruzese et al. (1984) and 
Strobel et al. (1985) to deduce an upper limit on D H of 6 x 10 5 cm s from 
thermodynamic considerations only. Model A in Fig. 4 is the appropriat 
height dependent upper limit on D H . A comparison of these results with 
earlier work of Johnson and Gottlieb (1970) based on the Kuhn and London 
(1969) IR cooling rates indicates agreement in the critical mesopause region 
to within 50%, the Johnson and Gottlieb D H values being larger. It should be 
noted that they did not include gravity wave heating in their calculation and 
other input quantities such as solar flux, IR cooling rates, 0 recombination 
heating rates, etc., were different from values used by Apruzese et al. (1984) 
and Strobel et al. (1985). But the Johnson and Gottlieb (1970) results are 
still an accurate representation of heat balance constraints on gravity wave 
induced diffusion. The gravity wave heating term (H, f Eq.2) is at most 25% of 
the total solar UV and O recombination heating rate in the calculations of 
Apruzese et al. (1984) and Strobel et al. (1985). 

To further illustrate aspects of the above discussion, the analytic 
expressions of Strobel et al. (1985) are adopted for the gravity wave terms in 

Eq. (2) with e - 1 
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( 6 ) 


",-c„-i.s(^*i-Aj Dv (, 0 < cmV ') unUs [ K d -,j 

where 

1 *D U 

SI ry 

D u dz 

When P r - 1 and D H is constant, then H g - C g - -7.5 and -22 K d' 1 for D„ - 10 6 
and 3 xlO 6 cm 2 s' 1 , respectively. These values should be compared to the 
total, globally averaged, solar heating rate including O recombination of 
Apruzese et al. (1984), 

z(km) — 65 70 75 80 85 90 95 100 

Quv(K d 1 ) - 3.2 1.6 1.7 5.0 9.2 8.9 13 18 

The upper mesosphere (70-80 km) is thus seen to be a critical region where the 
solar heating rate is low. With constant D„ (— D H ) a value as small as 2 x 10= 
cm 2 s' 1 is sufficient to balance solar heating at 75 km. In competition with 
solar heating is C0 2 infrared cooling primarily in the 15 n bands. The most 
accurate calculations of this cooling are by Dickinson (1984) who included as 

accurately as possible non-LTE effects. His globally averaged cooling rates 
are 

z(km) - 65 70 75 80 85 90 95 100 

Qir(K d 1 )-3.0 1.0 0.5 1.1 3.0 4.4 7.4 11.5 

A comparison of and below the mesopause suggests that gravity waves are 
not driving the globally averaged structure of the mesosphere far from 
radiative equilibrium. If one deduces values of from a balance of - Q 1R 
with H g - C g given by Eq. (6) one would infer mesospheric values similar to 
Model A in Fig. 4. These values up to 80 km would be consistent with the 
constraints on K II (— D H ) derived from ground-based microwave spectroscopic 
data discussed above. 
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One possibility for strong gravity activity in the mesosphere without 
substantial effects on the transport of constituents and heat is for the 
induced diffusion over macroscopic scale lengths to have a large effective 
Prandtl number, a point advocated by Strobel et al. (1985,1987) with strong 
theoretical support fro. Fritts and Dunkerton (1985). Fritts and Dunkerton 
(1985) examined constituent and heat fluxes driven by localized gravi y 
breaking where the breaking zones are small in vertical extent In comparison 
to the vertical wavelength. In their analysis localization of turbulence 
yields an eddy Prandtl number greater than 2 between saturation (A - 1 for 
saturation amplitude) and modest supersatur.tion (A - 1.3). Only for large 

supersaturation (A - 2) does P r approach 1. 

For constant D„ in Eq.(6), P, - 6 yields H, - C, - 0, i.e. no net heating 
or cooling from gravity waves. High effective eddy Prandtl number turbulence 
would allow significant deceleration of the mean zonal winds without a 
substantial diffusive signature in the thermal structure. The physical reason 
for this result is that the conversion of the gravity wave's kinetic energy 
still produces significant heating of the atmosphere although conversion of 
internal energy is reduced. High Prandtl number turbulence leads to sluggish 
eddy transport of heat and large reductions in the associated divergence of 
the eddy heat flux. As a consequence the conversion of gravity wave energy 
into heat can approximately balance the divergence of the eddy heat flux. 

It is also worthy to not. that Justus (1967) determined P r to be - 3 from 
photographic tracking of rocket released chemical clouds and analysis of 
turbulent wind data in the 90-110 km region. There was considerable scatter 
in his data points but according to his error bars P, was at least 2.2. 

Net gravity wave heating occurs also when D„ increases sufficiently 
rapidly with height (small H„) ; but at some altitude it must level off with a 
consequent large divergence in the eddy heat flux. The model results 
Apruzese et al. (1984) and Strobel et .1, (1985) yield gravity wave heating by 
this circumstance of. at most, 0.4 K d‘> and limited to below 80 km. 
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In Eq. (2) the term involving the divergence of < P wT> has been omitted in 
the above quantitative discussions. The only model results available to 
consistently ascertain its importance are from Garcia and Solomon (1985). For 
solstitial conditions at the mesopause with a mean circulation driven by 
gravity wave breaking, their globally averaged values are - D H - 1.5 x 10 6 
cm 2 s 1 , and <sr>-i6* e m *". The w and T fields have an approximate height 
dependence of ft. The convergence of this heat flux is spread over at least 
2-3 scales heights and has a magnitude of - 1 K d' 1 , substantially less than 
the value of H g - C g - 11 K d' 1 for this high diffusivity, P r - 1 model of the 
middle atmosphere. If the convergence of the < P wT> flux were this large for 
values of D H - 10 5 cm 2 s' 1 , then this downward heat flux would be have to be 
evaluated more carefully in globally averaged heat budgets. The fact that 

this heat flux is not large is not surprising as w and T have large amplitudes 
only at polar latitudes. 


5. Constraints on D H 

Of all the eddy diffusion coefficients discussed so far this is the most 
difficult to evaluate. The fundamental effect is the deceleration of the mean 
zonal and meridional winds. Theoretical models of the mean circulation of the 
middle atmosphere are sensitive to a number of input parameters. The 
radiative drive depends on the net imbalance of Q OT and Q IR . Whereas the 
former can be reasonably accurately calculated in a large dynamical model, the 
latter cannot be computed at the level of detail and accuracy of Dickinson's 
(1984) non-LTE model with present day computer resources. The infrared 
cooling acts also to damp temperature perturbations associated with the mean 
circulation and accurate simulation of the temperature and wind fields depends 
critically on the radiative damping. 

The deceleration of the mean zonal wind by gravity wave stresses is 
written as 


da 

— t- 

dt 


\ d . 

-- (PU „ ) 


A 

u - c 


(7) 
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Thus an approximately accurate value of D M can only be inferred from the 
deceleration if independent knowledge of the phase speeds of the gravity waves 
is available. Furthermore the right hand side of expression (7) holds only if 
a saturated gravity wave maintains constant amplitude with height, as Lindzen 
(1981) originally hypothesized. Recently this hypothesis has been questioned, 
Schoeberl (1988) argues that the amplitude of a saturated gravity wave can 
still grow with height. Only when the convection zone is comparable in 
vertical extent to the vertical wavelength will wave growth be significantly 
attenuated. 

On the basis of available estimates for D M from Lindzen (1981) , Holton 
(1982, 1983), and Garcia and Solomon (1985), it would appear that deceleration 
of the mesospheric zonal winds require that exceed 10 6 cm 2 s" 1 on a globally 
averaged basis in the upper mesosphere, a value definitely in excess of K zz 
and D h . For example, a typically deceleration rate of 100 m s' 1 d 1 (Holton, 
1983), would for (i-cl-wms' 1 yield — 1.3 x 10 6 cm 2 s 1 . 

6. Concluding Remarks 

The most stringent constraint on gravity wave induced diffusion is 
obtained from ground-based microwave spectroscopy measurements of H 2 0 with 
support from microwave measurements of CO and rocket-based mass spectrometer 
measurements of Ar (von Zahn, 1970). These data constrain K„ (- D H ) to be 
(1-3) x 10 5 cm 2 s' 1 in the upper mesosphere. Atomic oxygen data also limit K„ 
to - 10 5 cm 2 s' 1 at the mesopause, but indicate an order of magnitude increase 
in K,, z in the first 10 km of the lower thermosphere (Brasseur and Offermann, 
1986). 

The evidence cited above, which leads to the conclusion that K ZI is low 
in the mesosphere, is based on the absolute concentrations of the selected 
chemical tracers. An alternate point of view can be advanced, based on the 
seasonal variations of 0 and 0 3 . The green line radiated by 0 shows a 
semiannual variation which Garcia and Solomon (1985) argued can only be 
explained by a semiannual variation in the gravity wave induced diffusive 


12 



transport. The equinoctal periods are characterized by weak zonal winds and 
hence, by the Lindzen ( 1981 ) parameterization in their model, weak gravity 
wave induced diffusion. Soltistial periods have characteristically strong 
zonal winds and are accompanied by strong diffusion, thus creating the 
semiannual variation in diffusive transport. This variation in diffusive 
transport should also produce a semiannual component in the seasonal variation 
of the mesospheric H 2 0 mixing ratio profile. But Bevilacqua et al. (1987b) 
found from ground-based microwave measurements of water vapor no obvious 
semiannual component, only a pronounced annual component. Ozone mixing ratios 
inferred from SME data exhibit a pronounced semiannual variation particularly 
in the spring of the first two years of data acquisition (1982-1983) at 80 km. 
Data from later years do not contain such a distinctive component. The 
existence of this semiannual component may depend to some degree on vertical 
displacements of the ozone profile as it contains a minimum value in mixing 
ratio at — 80 km, which if displaced by a few kilometers could create the 
appearance of a time varying component at a fixed height with the periodicity 
of the vertical displacement. 

Consideration of the energy balance of the upper mesosphere indicates 
weak departures from radiative equilibrium, consistent with D H < 6 x 10 5 cm 2 
s" 1 at the mesopause and decreasing substantially with decreasing altitude. 

The uncertainty attached to the independent calculations of solar UV heating 
and atmospheric IR cooling cannot exclude the possibility of the upper 
mesosphere being in radiative equilibrium and the inferred values of D H being 
upper limits. Similarly the uncertainty cannot exclude a value of Dg as large 
as 10 6 cm 2 s" 1 . 

The available calculations to date suggest that D M is at least 10 6 cm 2 
s~ l . The consequences of an order of magnitude smaller momentum diffusion 
coefficient due to more modest reductions in zonal wind dfeceleration and in 
(u-c) have not been adequately explored, although we note that Schoeberl (1988) 
advocates reduced wave stress on the mean flow in the mesosphere. It is worth 
citing the lidar studies of the nighttime Na layer over Urbana, IL by Gardner 
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and Voelz (1986). They found average values for mean flow deceleration of 
-27.2 m s' 1 d _1 and a corresponding eddy diffusion coefficient, presumably D M , 
of 1.8 x 10 5 cm 2 s' 1 . Individual measurements yielded deceleration rates up to 
-200 m s' 1 d' 1 , but the long term averages were an order of magnitude less. 

But Reid and Vincent (1987) found the mean zonal flow deceleration to be 
typically between -50 and -80 m s' 1 , with occasional values as large as -190 m 
s" 1 . Also they measured comparable mean meridional deceleration rates. 

Likewise Meek et al. (1985) obtained comparable deceleration rates. Clearly a 
global climatology mean wind deceleration needs to be constructed from the 
available, but sparse data base and more precise altitude information is 
needed to determine where the deceleration actually occurs. For the present 
we are in somewhat of a dilemma; depending on what evidence is adopted, a case 
can be made for the Prandtl number (D„/ d h) applicable to turbulent mixing over 
macroscopic scale lengths to be anywhere between 1 and 10, although model 
values of D„ clearly suggest a large Prandtl number. In this author's 
judgment the available evidence suggests that the Prandtl number is > 3. 

Most of the inferences of eddy diffusion coefficients for constituent and 
heat transport were accomplished with ID models. The exclusion of mean 
circulation transport of constituents and heat may render the inferred values 
of gravity wave induced diffusion, in many instances, to be only upper limits. 
The low values deduced for K ZI and D„, in fact, suggest that the mesosphere is 
advectively controlled rather than diffusively controlled. 

Acknowledgments: This work was supported by NASA Grants NAGW-826 and 
NAG5-796 at The Johns Hopkins University. 
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8. Figure Captions 


Fig.l. Comparison of ground-based, observed H 2 0 mixing ratios with 1 a error 
bars (Bevilacqua et al . , 1985) and model results for the D H profiles 
given in Figure 2 and with chemical acceleration for 0 and 0 3 . After 
Strobel et al . (1987). Copyright by the American Geophysical Union. 

Fig. 2. Vertical eddy diffusion coefficient (D H ) profiles for indicated models 
in Figure 1. After Strobel et al. (1987). Copyright by the American 
Geophysical Union. 

Fig. 3. Calculated globally averaged temperature profiles for the D H profiles 
given in Figure 4, with P r — 1, and compared with the CIRA 1972 
temperature profile. After Strobel et al. (1985). Copyright by the 
American Geophysical Union. 

Fig. 4. Vertical eddy diffusion coefficient (D H ) profiles for indicated models 
in Figure 3. After Strobel et al. (1985). Copyright by the American 
Geophysical Union. 




(ran) apmniv 


Figure 1 


Mixing ra 





profiles 



Figure 2 


10000 IE-4-5 1E+6 1E + 7 

Eddy diffusion coefficient (cm a sec““ 1 ) 






TEMPERATURE 










DIFFUSION COEFFICIENT 



(UO|)Z 

Figure 4 


(i- Sjtuo) 


